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Abstract 

We present a calculation of the two-loop evolution kernels of the twist-two transversally polarized 
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and discuss a solution of the evolution equations suitable for numerical implementation. 
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1 Introduction. 



Helicity-flip parton densities, i.e. densities of transversally polarized quarks and linearly polarized 
gluons, are inaccessible in conventional deep inelastic scattering experiments on polarized spin-| 
targets. The transversally polarized quark density |I| is of odd chirality, and thus requires a helicity 
flip in a perturbative QCD subprocess which is forbidden for massless quarks, while the linearly 
polarized gluon requires the helicity to be flipped by two units which is not allowed by angular 
momentum conservation for spin-| hadrons. However, it will definitely appear for spin-J > 1 
hadrons/photon ^. 

If one goes from the forward to off-forward kinematics, i.e. inclusive DIS to exclusive DVCS 
with a real final state photon 0, |^, ^], one allows for non-zero orbital momentum in the system and 
thus the gluons can transfer two units of helicity from photons to spin-i hadrons [0, H, ^. Due to 



specific cos / sin (30) azimuthal angle dependences in the factorized p, |10|, [TT[] cross section of the 
corresponding skewed tensor gluon distribution (SPD), it can be extracted |^ without the usual 
complication due to quark contamination. A similar possibility exists for a system of hadrons 



due to their orbital motion, as it occurs in the 77 vrvr reaction As to the quark sector, 
unfortunately, the original proposal to measure the skewed quark transversity in diffractive meson 
production [^] has been proven erroneous due to the preservation of chiral symmetry for the cross 
section [T^l and thus complicates the experimental access to the function in question. 

Note, however, in DVCS, since the photon helicity flip requires a change in the hadron helicity 
by two units, the amplitude on a proton target is suppressed (due to its spin being only i) by 
an extra power of the transverse momentum transfer in the t-channel A_|_ whereas it will not be 
suppressed on a deuteron target, which is available at certain facilities like HERMES. Given its 
current run on high density targets, it is very interesting that high statistics data will be shortly 
forthcoming and hence there is a need to study the helicity-flip amplitudes theoretically. 

In this paper, we present an evaluation of the two-loop skewed evolution kernels for chiral 
odd quark and tensor gluon distributions. Since corresponding operators which define the non- 
perturbative content of the former belong to different representations of the Lorentz group the 
mixing is absent in both cases. Our presentations is structured as follows. In the second section, 
we discuss the general solution of the skewed evolution equations, in the third one, we present the 
reconstruction of the skewed evolution kernels in the Efremov-Radyushkin-Brodsky-Lepage (ER- 
BL) and Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) regions, and finally, we summarize. 
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2 Evolution equation. 



SPDs IQ-P, appear in various exclusive processes and their extraction would provide us 



new, valuable, non-perturbative information on the structure of hadrons. These functions can 
be considered as a link between Feynman's parton densities and exclusive quantities like form 
factors/distribution amplitudes, in fact they are generalizations of both. In QCD they are defined 
as expectation values of twist-two light-ray operators over the hadronic states, 

{ ) Q') = { 4P-. ) / { } --^'i^-e'l^'). (1) 

where P+ = n - (Pi + P2) with being a light cone vector. These are functions of the momentum 
fraction x, the longitudinal momentum fraction 77 > in the t channel, called the skewedness 
parameter, the momentum transfer square = (P2 — -Pi)^; and they vary depending on the 
resolution scale Q^. The partonic interpretation of these functions depends on the kinematical 
region. For the ER-BL region |x| < r/ they behave as distribution amplitudes of meson-like states, 
while for the DGLAP domain |x| > 1] they mimic the conventional probability distribution of 
finding a parton with a corresponding momentum and spin state that is scattered from the initial 
to the final state. 

Now we recall the general properties of the evolution equation for SPDs, 
d 



(2) 



d\nQ^ 

with the kernel ^^V^ defined as a series in the coupling as- Note, that we treat both quark 
and gluon cases simultaneously and that there is no mixing between them, as mentioned in the 
introduction, since they belong to different representations of the Lorentz group. 

The evolution kernel ^^V^{x,y) with < x,?/ < 1 is of the ER-BL-type and governs the 
evolution of meson distribution amplitudes (formally \ri\ = 1). The general structure of this kernel 
was already studied in the past |Q with the main result being that the evolution kernels of skewed 
parton distributions can be obtained from the ER-BL ones, 

^^V^ = e{y - x)f{x, y) + e{y - x)g{x, y) + T ^ "l] , x = l-x, < x, y < 1, (3) 

{y ^y] 

[here x denotes a new variable, which differs from that one used in Eq. (^] by replacing the 9 
structure: 

for X — > ^(r] + x), and analytical continuation of the prefactors / and g. In the gluonic case 
the second structure can be reduced to the first one by means of Bose symmetry, namely. 
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'^q^{x,r]) = '-^q^{—x,r]), while in the quark case, we should decompose the SPD in contributions^ 
coming from quarks {x > 0) and anti-quarks (a; < 0). 

The evolution equation (|]) can be solved by means of different methods. However, most of 
them are plagued by either a limitation of their applicability to a leading order (LO) analysis or a 
lack of sufficient accuracy in handling particular regions of the phase space. Only the orthogonal 



polynomial reconstruction method |T5|, |I6[ and the direct numerical integration methods [18 
allow for a generalization beyond the one-loop approximation. However, the first method suffers 
from numerical complications in the treatment of the transition region |x| = and, therefore, one 
may explore the efficiency of a direct numerical integration to overcome this particular difficulty. 
To do so we should separate the evolution equation into DGLAP and ER-BL regions. 

In the quark sector, it proves convenient to introduce (non-) singlet like combinations of SPDs 
that are symmetric/antisymmetric w.r.t. the momentum fraction x: 

%'^{x,ri) = ^ ^q^{x,ri) — ~q^{—x,ri) , '^q'^{x,ri) = ^q^{x,ri) ^%'^{—x,ri). (5) 

Here the antiquark contribution is given by —'^q'^{~x,ri) and thus ^q'^{x,ri) for x > is just 
the sum and difference of quark and anti-quark distributions, respectively. The kernels in the 
corresponding evolution equations read 

QQ r± ( v + x v + y \ v + x v + y \ qqpt± ( v + x v + y \ jx^~x\ 

[ 2v ' 2v J 2v ' 2v J [ 2r^ ' 2v )^\y-^-yj' 

Based on the support and symmetry properties, we discuss now the solution of this evolu- 
tion equation (^. In the following ?7,Q^) is a (anti-) symmetric skewed distribution of any 
species, where for gluons we formally have to take the solution for ~q{x, rj, Q"^). An analogous sign 
convention holds true for unpolarized partons, while for longitudinal polarized ones, one should 
replace + ^ — . In the DGLAP region we have a homogeneous integro-differential equation, 

while in the ER-BL region an inhomogeneous term arises, 
d ± ( .r^'i^ _ /■''^' dy ±.Jr] + x r] + y 



dy , r] + x r] + y\^ , ±Tf 



±F ( ^ + ^ v + y \ -^±p (v-x r] + y 



21] ^ 21] J \ 21] ^ 21] ^ 

^More precisely, we should speak of (anti) quark distributions only for x > rj and (a; < —77), while there are two 
meson-like distributions with different symmetry in the ER-BL region. 
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where we used the short-hand notation F{x,y) = F{x,y). Introducing the kernel P (^|, = 
(^F — (^2^' ^2^)' obvious that equation (|^) is of the DGLAP type. Since, in addition, 
the kernel depends on the skewedness parameter t], there does not exist the same factorization for 
the Mellin moments as in the forward case. Fortunately, this equation can be efficiently integrated 
by brute force . We formally write the solution of this equation as a convolution of the skewed 
distribution given at the input scale Ql with an evolution operator 

g(a;,?7;Q^)l^.|>^ = / —UY)GLAp[x,y,v;Q^,Ql)q{y,V,Ql), (9) 
J\ri\ y 

that satisfies the equation 

d dy fx ^\ 

-T{—^^^GLAp{x,z,r];Q^,Ql) = / —P[-,-]UBGLAp{y,z,r];Q^,Ql) (10) 
amc^ Jx y \y y J 

with the initial condition [/dglap (a;, y, r]; Qq) '^o) — '^(-'-^f)- '^^^ homogeneous part of the 
evolution equation (H) is precisely of the ER-BL type and its solution is 

g^°"^(x, Q2)i.|<, = £ dy L/erbl 2/; Q\ Ql^ g'°'"(w, r/; Ql), (11) 

with U-EKBh{x,y,Ql,Ql) = 6{x — y). A particular solution for the inhomogeneous equation is 
easily constructed by means of the inverse evolution operator and we obtain, together with the 
homogeneous part. 



g(x,r7;Q^)l^l<^ = J^^dy Uerbl (^^,y;Q^,Ql^ ^q{vy,V;Qo) 

+ j ^ dz j^^ f/Ei^BL (y, Q'\ Q^) /(t^z, r/; g'2)|, (12) 



where /(x, r^; Q'"^) is the convolution of Eq. (|^) with the solution of Eq. (^) in the DGLAP region. 
An important issue is to study the evolution at the point \x\ = ri. Since we have the property 



F{x = 0,y) = F{x = 1,?/) = in both LO and next-to-leading (NLO) [see below Eqs. (|T8 



and (^01)], the limit e ^ of the evolution equations for x = r] + e and x = r/ — e will be the 
same. Therefore, one expects that a smooth input function will remain smooth under evolution. 
Nevertheless, this point necessitates care as it is a delicate point in the numerical solution of the 
evolution equations. 

The construction of the ER-BL from the DGLAP kernels is in principle possible, provided we 
know the eigenfunctions of the former. In LO the exclusive kernels are diagonal w.r.t. Gegenbauer 
polynomials 

dx ^(A) ( x\ AA^mT ( v + ^ v + y \ _ ^AAjo)T^u(A) ( y\ , . 
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where iy{A) = {§, §} for A = {Q, G} and ^^f^ are the forward anomafous dimensions at LO. 
This fact reflects the underlying conformal symmetry at tree-level. Consequently, this symmetry 
allows us to reconstruct the ER-BL kernel from the DGLAP kernel by an integral transformation 
[for details see second article of [g], ^^V^^^^{x,y) = dz G{2x - 1,2?/ - 1; z\u{A))^^p(^^^{z), 
where the integral kernel is given in terms of a hypergeometric function 

i (-)l(2z/) [l - 2 [xy - ^{1 - x2)(l - y^)) z + z^ 



•-2' 



1^ + 1,1^ 

2u 



A^{l-x^){l-y^)z 



\ 



l-2(xy- ^(1 -x2)(l -?/2)j z + z'^ ^ 



Note that, in any order of perturbation theory, the evolution operator in the ER-BL region can 
be represented in terms of an infinite series of Gegenbauer polynomials [^. Unfortunately, since 
in general a SPD does not vanish at the point x = ±7], it is expected that, to achieve a good 
approximation, a large number of terms in this expansion is needed. Fortunately, we can resum 
the Gegenbauer expansion in the ER-BL region which provides us the integral transformation, 
mentioned above, of the evolution operator in the forward region 

1 

f^ERBL (a;, Q^ Qo) = 2 Jo '^^ ^' ^I'^'^'^^^^dglap (2;, 2;' = 1, = 0; (5^ Qq) , (15) 

where G{x,y; z\v) is the kernel defined in Eq. ([1^ . Beyond the LO approximation, one can per- 
turbatively expand the solution with the help of this evolution operator. From the representation 
of this operator in terms of Gegenbauer polynomials it is obvious that the inverse operator is 

f^ERBL (a;, y] Q^ Ql) = %rbl (x, y, Q^). 

Let us add that in the forward case, i.e. (P2I — Eqs. (|l|) define the chiral odd quark 
density and the tensor gluon density times the momentum fraction z, thus, we have the limit of 
the splitting functions: 

I GGpT I (^) = LIM I ^^^^ I (x, y) = lim — | j ( ^ ^ j ' ^^^^ 

which will be of relevance to the analysis in the following section. 



3 Construction of two-loop kernels. 

Since the extension of the ER-BL kernel to the whole region is unique , it is sufficient to construct 
the evolution kernels, 

^^V^{x,y) = ^^^V^^'\x,y) + (g)' ^^K^W(a:, y) + O (a^) , (17) 
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in the ER-BL-region. The one loop kernels have been obtained by a direct calculation pll, |22 



Alternatively, they can be deduced by means of conformal covariance ([T^) from the forward kernels 



via the integral transformation (14): 



( QQyTiO) ^ 
j GGyTiO) 

where 

QQjT ^ QQj!b ^ 



e{y - x) 



T • 



GGj^T 



X 

y 



X 

y 



GG^no) 



6{x - y), 



X 



yy 



X 



GGrT 



r 



GGrb 



r 



X 



y^ y 



X 



GGTiO) 

7i 



6Ca + Po. (19) 



Here the +-prescription is defined for both channels in the same way, namely as [V(x, ?/)]-|_ = 
V{x, y) — S{x — y) j\ dzV{z, y). As explained in great detail in in NLO conformal covariance 
does not hold true anymore for the minimal subtraction (MS) scheme in the dimensional regular- 
ized theory. Besides a symmetry breaking term proportional to the QCD /3-function, there appears 
an additional term which is induced by the leading order anomaly in the special conformal trans- 
formations of conformal operators. Note, that the latter anomalies are renormalization scheme 
dependent even in one- loop approximation. Next, the results for the local conformal operators in 
the MS scheme can be transformed to the kernels, which have the following structure in NLO: 



AA 



AAyT 



AAyT(Q) ^ j _ 



AAgT^ 



AAyT{0) 



AA^T 



(20) 



where the commutator stands for [A0 B]{x,y) = Jq dz {A{x, z)B{z,y) — B{x, z)A{z,y)}. The 

off-diagonal part w.r.t. Gegenbauer polynomials is contained in the first two convolutions on the 
r.h.s. of this equation. The third term, i.e. "^^V^, is diagonal w.r.t. Gegenbauer poljTiomials 
with index z/(y4) = ||, || for A = {Q,G}. The dotted kernel is derived from the LO ones via a 
logarithmic modification. 



QQyTiO) 
GGyTiO) 



Cpdiy 



X 



XI 



CA0{y-x, 



yy 

x"^ 



X y 



X 

y 



X 

y 



In - + 



X 

y 



y y -X y 

The (7-kernels, meanwhile, contain new information and are defined by 



(21) 



f22) 



QQgT 



GGgT 



-Ca 



e{y 



e{y 



X 



XI 



In ( 1 - ^ 

V y 

y - X 



+ 



y-x 



X 

[y 

- 2- 



X 

y 



(23) 
(24) 
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In contrast to the QQ-kernel, which is unique in all sectors |]T9[, the GG-kernel differs from those 
in the chiral even case by the term 2C^^ 0. 

Since it is difficult to project onto the off-diagonal part and apply the integral transformation 
to the NLO DGLAP kernel we use another strategy for the construction of the remaining 
diagonal term, T>. This term is decomposed into a contribution related to the crossed-ladder 
diagram containing dilogarithms, G, and a remainder D which is, for the present helicity-fiip 
sector, given as a linear combination of one-loop kernels. 



jQQV^\ 1 



CpiC. 



In the following the kernels are defined as 



r QQjjT ^ 
[ GGDT ] 



+ <! ™^ ) {x,y). 



(25) 



AA^T ^ ^AA^T\^ y)j^e{y~x) (^^F + A^^r) (x, y) , (26) 



where 



AAj^T 
AA-^T 



2^^/^lnxln?/-2^Y^[Li2(x) + Li2 



AArT 



(27) 



AAfT 



r 



2Li2 1 



X 



- y 

y 



+ 2 [Li2(y) - In X In y] + 2 """"fU^^x 



AAfTi 



Here the QQ kernel follows from the explicit flavor non-singlet two-loop result in the chiral even 
sector |2|] by replacing QQ/^s _^ QQjt ^ ^j^j^ addenda A^^h'^ and A'^'^h'^ being zero. Note that 
we rewrote the original result in such a way that both ^-contributions are symmetric with respect 
to the weight function xx, i.e. yy^^h'^{x,y) = xx^^h'^{y,x) and yy^^W{x,y) = xx'^'^W{y,x). 
Now it can be verified that both parts are separately diagonal w.r.t. Gegenbauer polynomials. 

Since the two-loop crossed-ladder diagram contains only a simple pole in the dimensional 
regularization parameter, the A/" = 1 supersymmetric constraints on the level of local operators 
lead to the equality of chiral odd quark and tensor gluon anomalous dimensions '^^jj'^ = ^'^'jj 
From this constraint one can derive a relation between the momentum fraction kernels 



21, 25 



which reads [for details see p6[| ]: 

d 



d 



^ ^<^G^{x,y) + ^^^G^{x,y) 
oy ox 



0. 



(28) 



Taking the generic form ( PBD for ^^G^ {x,y), this differential equation and the symmetry require- 
ments imply the following form of the gluonic addenda (the remaining degree of freedom is fixed 
by the diagonality of the lowest moments): 



A^^/.^(x,y) 



2x 

y'^y 



2x , _ 2a; , 
Inx — Iny, 

y y yy 



A^^F(x,y) = A^^'/i^(x,i/) 



GGuT/ 



(29) 
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The missing diagonal "^^D^ terms are easily obtained in the DGLAP representation by taking the 
forward limit and comparison with the two-loop results namely, 



-G-terms 



(30) 

Since the result is expressed in terms of LO splitting functions, we can easily restore the ER-BL 
representation from themQ: 



QQ„.a. 



(31) 



GG„.a I GG^c 



(X, y) + PoC^ { I'^'^v^ - I'^'^v' + \^^v^ \ (X, y) 



+ d \ - + ^^S'^ \ (a;, y). 



where the 6-kernel has been defined above in Eq. ([T^) and the a- and c-kernels are given by 



GG^a 



X 



GG^c 



X 



— {2xy + y -x) . 



(32) 



y y y 

Contributions concentrated in x = y will be restored below. 

It remains to perform the exclusive convolution. For two non-regularized kernels given by their 
functions fi{x, y) with i = {1, 2} we have 



{vi ® V2) {x, y) = 9{y - x) (/i (g) /a) {x, y) + 



y^y. 



(33) 



where 



ry /■! rx 

{fi® f2){x,y) = dz fi{x,z)f2iz,y) + dz fi{x,z)f2{z,y) + dz fi{x, z)f2iz,y). (34) 

Jx Jy Jo 

The convolution of two regularized kernels is slightly more involved. Here we define: 



K]+ ® N + ) {x,y) 



d{y~x) (/i®/2) {x,y) + 



where the convolution on the r.h.s. can be decomposed into well-defined integrals 



\ X - 


X 




\v- 







(35) 



(/i®/2) {x.y) 



dz{[fi{x,z) - fi{x,y)] 1/2(2, y) - f2{x,y)] + [fi{x,z) - fi{z,x)] f2{x,y)} 



+ / dz{[fi{x,z) - fi{x,y)]f2{z,y) -[fi{z,x) - fi{x,y)]f2{x,y)] 



+ / dz{[fi{x,z) - fi{x,y)] [f2{z,y) - f2{x,y)] + [fi{x,z) - fi{z,x)] f2{x,y)} 
Jo 

- fi{x,y)f2{x,y). (36) 



■^Note that the last term in the QQ channel was missed in the second article of Ref. 
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Performing the exclusive convolution in Eq. ( |20D and adding the diagonal terms given in Eqs. 
(^51) - (^71), (^), and (|3l|), the QQ and GG transversity kernels in NLO are found to be: 



QQy{l)T± 
GGy(l)T 



where 



Inxlnx 

y 



7r"5(a;-y), (37) 



I yy 



m X H — m X m a; 



hi 


X - 


X 






[y- 







i^^7r5(x-y), (38) 



cl< 



-GfIG 



- C(2)) 



QQfT QQrT 



Ca 
2 



'3 
^4, 



In — I In — 



/^)ln-ln(l- 



X 



y , 



'1 ip x^QQj" _|_ _QQjT _|_ 
3 



x,y)T^^h^{x,y) 



^^^^(x,l/) 

= Cl|2 (1 - C(2)) ^^/^ + ^^f In^ ^ - (^^/^ - ^^F) In ^ In (^1 - 



+- 
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3 2 



i^^/^^(x,y)-l^^F(a:,y)-A^^/.^(x,y) 



^GGjc 



I ( 



2 ' " 1^4 " 3 " 2 " J ^ / ^1 3 ^ 



X 



(39) 



(40) 



The contributions concentrated on the diagonal x = y are fixed by the lowest conformal moment. 
We have in the QQ and GG channels: 



GGn)T 
71 



12 ^ 12 



^ + (l±l)[7-8C(2)+4C(3)] 



(41) 



Finally, we derive the NLO skewed DGLAP kernels in Radyushkin's notation with fraction 
z = (here x denotes the momentum fraction in the SPDs) and skewedness ( = These 
kernels govern the evolution of SPDs in the kinematic range of C < z < 1 and —l + (^<z<0. 
We choose this set of variables such that we have the closest resemblance to the usual DGLAP 
kinematics. In fact for DVCS and vector meson production ( = xbj up to terms ~ ^! The 
derivation procedure follows that of the forward limit reduction of the extended ER-BL kernels. 
After the analytical continuation of the ER-BL kernels to the whole {x, y} plane by the replacement 
of the ^-function structure, we derive the kernels in the following way which is suggested by the 
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non-zero support of the ^-functions after the replacement >1 — >^ and y — > 1 — ^ 

for the above mentioned kinematical regime: 



c 



AAyT 



AAyT 



(42) 



Note again that we omitted the additional factor of z^^ in the GG kernel so as to take into account 
that in the forward limit, the skewed gluon distribution turns into Q^) rather than Gi^z^ Q"^)- 

The results for the chiral odd DGLAP kernels in LO, following our above prescription, are: 



QQpT^ 



1 + z 1 



GGpT ^ {^ + Z)Z 1 + C + ^ 



Z 



e c 



(43) 



and in NLO: 



1 



ln(l-gln 



+ - In 2; I In ^ 

z 



+ In 2; In ^ 



1 1 



+ ( - - 2C(2) ) 



- 7;CfPo{ Ti'^'^p' +^lnz+ 1 - 9 In 1 - 9 



Cp c 



(1 T l)Wp« + ^QQpT + QQhT ^ qq-^t]^^ 



(44) 



GG 



'cGpT _ ^ + C(_1 + ^) ^ + 4^ + 



c 



"^'^ f }_GG c I ^GG a 

4 V3 ^ 2 ^ 



2$ f C7 + 7 In ^ - C fl - 5 fl - 1 + - ^ + ( ^ - 2C(2) ) ^^p^ 



cv C c 

fin f 1 



c 



3 



GG^T 



^, +ln5lnzj -^^^/i^-^«^F 



GG^a I ^GG^c 



(45) 



where one has to regularize the whole kernel with the +-prescription and add the end-point 
contribution —^^^q^'^'^ 5{1 — z) and — ^'^'-^7^ ''"^5(1 — z) specified earlier. The ^"^h^ and "^"^h^ 
functions are given by the equations 



^^/i^ = 2^V(^,C) 
+ 2^^^f(^,C) 



InCln ( 1 - ^ 1 + Li2 (l--\-U,(l--^\- C(2) 



InC 



2 Lio 1- - -Lio 1- - 



AATT 



r = 2^V(C-^,C) 



InCln 1 



1 



InC In 



-|- In 2; — In 



Lis 1 



-Li2 ( 1 - 7 ) - C(2) - ln(C + ^) In 1 - f + In z - Li2 (C - ^) - Li2 -I 
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In 1 



In 



InC -2 Li2 1 



Li, 1 



The remaining skewed kernels appearing in the above formulae read 



QQpO. 



GG^a 



2— - ^(l- z'^) 



GG 



P =T-2 



and the "^"^kj were found to be: 



(46) 



(47) 



Z 



Note, that we have included the addenda of the skewed G-f unctions in the single logs and the 
rational functions of the respective colour structures in the kernels so as to simplify the overall 
presentation of the results. 



GGuT 



z 



GGuT 



(4J 



4 Conclusion. 

In conclusion, we have presented a calculation of the two-loop exclusive evolution kernels for 
skewed chiral odd quark and tensor gluon distributions, thus allowing for NLO analysis of the 
cross sections where these functions contribute. Our modus operandi is based on the formalism 



developed previously in Refs. |^6[ for vector and axial sectors. Due to chirality conservation only 
the helicity flip gluon SPD shows up in the DVCS cross section. A peculiar azimuthal angle 
dependence of asymmetries on the nucleon target will give an opportunity to pin down the 
magnitude of the gluon SPD provided high accuracy data will be available. This will shed some 
light on the unknown tensor gluon content of the nucleon or higher spin hadrons. 
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